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Abstract 

We consider a simplified model for Majorana fermion dark matter and explore constraints 
from direct, indirect and LHC collider searches. The dark matter is assumed to couple to 
the Standard Model through a vector mediator with axial-vector interactions. We provide 
detailed analyses of LHC mono-jet searches and IceCube limits on dark matter annihilation in 
the Sun. In particular, we develop a method for calculating limits on simplified WIMP dark 
matter models from public IceCube data, which are only available for a limited number of 
Standard Model final states. We demonstrate that LHC and IceCube searches for Majorana 
dark matter are complementary and derive new limits on the dark matter and mediator 
masses, including in addition constraints from LHC di-jet searches, direct detection and the 
dark matter relic density. 
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1 Introduction 

Weakly interacting massive particles (WIMPs) are attractive candidates for dark matter and 
are predicted by various extensions of the Standard Model (SM). Searches for WIMPs at the 
Large Hadron Collider (LHC) and through direct and indirect detection experiments are com¬ 
plementary, and probe different types of dark matter models and different regions of the model 
parameter space. 

To explore the nature of dark matter, and to be able to combine results from direct, indirect 
and collider searches, one may follow a more model-independent approach. So-called simplified 
models (see e.g. [1-3] and references therein) describe dark matter and its experimental signa¬ 
tures with a minimal amount of new particles, interactions and model parameters. They thus 
allow us to explore the landscape of dark matter theories, and serve as a mediator between the 
experimental searches and more complete theories of dark matter. 

Many simplified models for dark matter have been proposed in the literature. Minimal mo¬ 
dels describe dark matter by a single particle which interacts with the SM through a single 
mediator. We focus on a model with Majorana fermion dark matter and a vector mediator with 
axial-vector couplings to quarks. Such models predict spin-dependent WIMP-nucleon scattering 
cross sections and can thus be probed by IceCube in the search for dark matter annihilation in 
the Sun [4, 5]. Previous work on dark matter models with (axial-)vector mediators [6-10] has 
mainly focused on Dirac fermion dark matter and, in particular, has not analyzed dark matter 
limits from IceCube. 

The complementarity of IceCube and LHC searches for dark matter has been explored in 
Ref. [4], albeit in an effective field theory (EFT) approach, where the masses of all particles 
mediating the interaction between the SM and dark matter are assumed to be large compared 
to the energy scale of the process. The EFT limit is reliable to interpret the low-energy WIMP- 
nucleon interactions which are probed by direct detection experiments, but it may break down 
when analyzing dark matter searches with IceCube and with the LHC. While the capture of 
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WIMPs in the Sun is well described by an EFT, the annihilation process is in general not. 
Indeed, as we will show, the IceCube limits are sensitive to details of the simplified model which 
cannot be described within the EFT. Furthermore, the EFT may break down when probing 
dark matter production at the LHC: when the energy scale of the interaction is near or larger 
than the mass of the mediator, resonance effects become important, and the mediator has to be 
included in the particle spectrum of the model [11-20]. 

We thus provide a comprehensive analysis of collider, direct and indirect detection constraints 
on Majorana dark matter with vector mediators in the simplified model framework. The paper 
is organized as follows. The dark matter model is introduced in section 2. We derive limits on 
the model parameters from LHC searches in section 3. In particular, we use the most recent 
results for mono-jet searches from both the ATLAS [21] and CMS [22] collaborations. Section 4 
addresses constraints from the dark matter relic density, from direct dark matter searches, and, 
in particular, searches for dark matter annihilation in the Sun with IceCube. For a mediator 
that is lighter than the WIMP, annihilation into a pair of mediators can be dominant. In order 
to determine the IceCube model rejection factor in this region of parameter space we develop 
a method to estimate the limits for annihilation into two mediators on the basis of the limits 
for annihilation into the SM particles the mediator decays into. We validate this method by 
applying it to annihilation into top quark pairs for which we find very good agreement with the 
most recent public IceCube dark matter annihilation limits [23] for that channel. We finally 
combine the various constraints and derive new limits on the dark matter and mediator masses 
in section 5. Our conclusions are presented in section 6. 


2 A simplified model for Majorana fermion dark matter 


We focus on a minimal model with Majorana fermion dark matter, X) and a vector mediator, 
Vji, with axial-vector couplings to quarks, 

£ D 5^ + 9q ■ (i) 


We assume universal couplings of the mediator to the SM quarks and neglect couplings to 
leptons.^ Thus the model has only four independent parameters; the couplings of the mediator 
to the dark matter and the SM quarks, g^, and Qq, and the dark matter and mediator masses, 
and My, respectively. We require both couplings to be gx-,9q < dvr. 

The axial-vector interaction of Majorana dark matter with nuclei leads to spin-dependent 
scattering cross sections and contributes significantly to the dark matter capture rate in the Sun 
(see section 4.2). Searching for dark matter annihilation in the Sun with the neutrino telescope 
IceCube can thus place strong limits on such models, which are competitive with direct detection 
bounds and with dark matter searches at the LHC. 

The width of the mediator, Fy, is determined by the particle masses and the couplings: 
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where niq- is the mass of the SM quarks. 

The dark matter cross section at the LHC and the WIMP-nucleon scattering cross section 
depend on the product of the couplings and on the width of the mediator, ^J9q9x and Ty, 
respectively. We will thus present our results in terms of Fy and y/gq9x: rather than in terms 
of the individual couplings gq and gx- 

^The case of different couplings to up- and down-type quarks has been considered e.g. in Ref. [10]. However, 
as discussed in Ref. [24], gauge-invariance sets very tight constraints on the difference of these couplings. 
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Figure 1: Non-accessible regions in the Mv/mx-Fv/Mv plane for x/9x9q ~ ^ shaded region), 0.2 (blue shaded 
region) and 0.1 (green shaded region). The width of the bands marking the boundaries of the shaded regions 
indicates the explicit dependence on nix'- the upper and lower edge of the band correspond to the boundaries for 
rux = 1 TeV and 1 GeV, respectively. 


The relation between Fy, x/9q9x individual couplings gq and provides some insight 

into the phenomenology of this model. In the region My > 2m^, and for any given value of 
the product of the couplings, x/gq^, we encounter a minimal value for the mediator width, 
below which there is no solution for the individual couplings within the model. These regions 
are shown in Fig. 1 in the My/m-^yFy/My plane. The unaccessible regions only depend very 
mildly on m^, as indicated by the width of the bands marking the boundaries of the shaded 
regions (here was varied between 1 GeV and 1 TeV). In the allowed part of the region where 
My > 2m^ there exist two solutions for gq for any given xj9q9x- To derive conservative limits 
on the model from di-jet production (see Sec. 3.2), we adopt the smaller value for gq in our 
analysis, unless this would cause g^ > dvr for a given ^J~g^g^. 

In the EFT limit My 3> m^, where dark matter and SM quarks interact through a 4-fermion 
operator with coefficient 1/M^ = gqg^/My, we find that Fy/My > 0.3 g^gq = 0.3 (My/M*)^. 
The LHC Run I data probe suppression scales M* < 1 TeV, see Sec. 3. Thus, in the region 
My 3> x/s where the EFT is valid, Fy/My is typically larger than one, inconsistent with a 
particle-like interpretation of the mediator (c/. Ref. [6]). 

3 Collider limits 

3.1 Mono-jet limits 

Weakly interacting dark matter particles can be detected at the LHC through their associated 
production with jets, electroweak bosons or heavy quarks. The search for such signatures to¬ 
gether with large missing transverse energy (MET) has been performed at the LHC Run I and 
is one of the central goals of LHC Run H [3]. In the following we will focus on signatures with 
mono-jets and MET as presented in [21, 22]. Searches for electroweak gauge bosons with large 
MET are important in general, but provide weaker limits for the dark matter model we consider, 
see e.g. [25, 26]. 

To simulate the experimental signature for our model, we have generated events using Eeyn- 
Rules 2.1 [27], Madgraph5_aMC@NLO [28] and PYTHIA 6 [29], including QCD processes 
with one and two jets in the hard scattering. We have used DELPHES [30] for the detec¬ 
tor simulation and implemented the cuts employed in the ATLAS [21] and CMS [22] mono-jet 
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searches in an in-house program. With the observed and expected number of events provided 
by Refs. [21, 22] we are thus able to set exclusion limits at 95% confidence level (CL) on the 
different parameters of the model. 


^/9x9q = 1, Cy = O.OlMv 


=1, Tv = 0.5Mv 



x/9x9q = 0-2, Cy = O.OlMy 


x/gx9q = 0-2, Cy = 0.5Mv 




Figure 2: Lower exclusion limits in the m^-Mv plane at 95% CL for the ATLAS (blue lines) and CMS (red 
lines) mono-jet searches. The limits for the simplified model (solid lines), for the EFT (dashed lines) and for 
the EFT applying the Q-truncation (dotted lines) are shown. Four slices of the parameter space: ^jQxQq = 1 , 
Fv — O.OlMv (upper left panel), yJgxQq = L W = 0.5Mv (upper right panel), x/9x9q = 0.2, Fy = Q.QlMy 
(lower left panel) and x/QxQq = 0.2, Fy = 0.5My (lower right panel) are displayed. The blue shaded region in the 
left panels represent the parameters space not allowing a consistent solution for the mediator width as a function 
of My,mx,x/g^- 


We have also studied our model in the EFT limit, where the interaction is described by a 
higher-dimensional operator of the form {g^gq/My) xlgl^X As the EFT is valid only for 

energy scales below the mediator mass, it has been proposed to restrict the momentum transfer 
in the s-channel, Q < My, when calculating cross sections at the LHC [12]. 

In Fig. 2 we show exclusion limits on the dark matter and mediator masses, for scenarios 
with a small or a large mediator width, Fy/Mv = 0.01 and 0.5, and small or large mediator 
couplings, x/gx9q = 0-2 and 1, respectively. Note that we consider Fy/M = 0.5 a rather extreme 
benchmark case, which however is commonly adopted in the literature. For even larger Fy/M 
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the narrow-width approximation to the cross section calculation may not be reliable, and the 
interpretation of the mediator as a resonance becomes doubtful. The results have been obtained 
for the ATLAS [21] and CMS [22] mono-jet searches interpreted in terms of the simplified 
model, the EFT, and the EFT with a truncation Q < My- Note that not all combinations of 
the parameters m^, My, Ty and y/g^Jjq are viable, as discussed in section 2: for a small width, 
Fy/My = 0.01, most of the region My > 2m^ is theoretically inconsistent. On the other hand, 
for a large width Fy/My = 0.5 the whole parameter region is allowed. 


CMS mono-jet, y/gx9q = 0.2, rrix = 15 GeV 



Figure 3: Theoretical prediction for the total production cross section (solid lines) and its observed upper limits 
(dashed lines) in the EFT limit with (green lines) and without (black lines) Q-truncation as well as in the simplified 
model with IV = O.OlMv (red lines) and IV = 0-5Mv (blue lines). Example for y/FJFi = 0.2 and rUx = 15 GeV 
for the CMS analysis. 

We hnd large differences between the EFT and simplified model limits. For My < 2my- the 
EFT limit (without Q-truncation) extents to much larger values of and drastically over¬ 
estimates the sensitivity in this region. By construction, this is not the case for the EFT limit 
employing the truncation on Q. Furthermore, also for My > 2m^ we find significant differences. 
This is due to the fact that My lies in the range of accessible energies at the LHC. The cross 
section can be greatly enhanced due to contributions from on-shell mediators. This leads to a 
larger sensitivity in the simplified model than in the EFT and is most pronounced for a small 
mediator width. However, for small My - smaller that the minimal MET of approximately 
150 GeV required in the search - the cross section in the simplified model decreases again and 
the sensitivity becomes weaker. This effect can be seen in the slice of parameter space with 
y/gx9q = 0-2 and Fy = 0.5My (lower right panel), where we cannot constrain the dark matter 
and mediator masses from the CMS analyses interpreted within the simplified model while within 
the EFT (without Q-truncation) a limit on My around 200 GeV is obtained. 

To illustrate the origin of the differences we present the cross-section predictions together 
with the observed upper CMS limits for the simplified model, the EFT and the EFT with 
truncation in Fig. 3. The results are shown exemplarily for x/g^Qq = 0.2 and = 15 GeV, and 
for small and large widths, Fy/My = 0.01 and 0.5. The observed upper limits, which involve 
the observed number of events, the background expectation and the signal efficiencies, do not 
depend significantly on the mediator width, and are similar for the simplified model and the 
EFT with truncation. However, the cross section predictions from the simplified model and the 
EFT differ by orders of magnitude, in particular for small mediator widths. Furthermore, the 
simplihed model cross section for ^/g/JF/ = 0-2 and Fy/My = 0.5 never exceeds the cross section 
limit, so that no constraint on the dark matter and mediator masses can be derived. Note that 
the slopes of the cross section limits and the cross section predictions are quite similar, so that 
small changes in the observed number of events, or small statistical fluctuations in the estimate 
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of the efficiencies and cross sections, may have a visible impact on the limits. 

We conclude that for our dark matter model, the EFT interpretation of the experimental 
LHC searches is only reliable if the mediator mass is larger than the accessible LHC partonic 
energies. However, as the current LHC limit on the suppression scale M* = My / yj9x9q is of 
0(1 TeV), an EFT interpretation would require extremely large couplings. On the other hand, 
the EFT limit with a truncation on Q results in overly conservative limits which do not fully 
exploit the potential of the LHC searches. We will thus use the limits obtained in the simplified 
model in the compilation of bounds presented in section 5. 

3.2 Di-jet limits 

Limits on our simplified model can also be derived from searches for the mediator particle in 
di-jet signatures. A detailed and comprehensive study of di-jet constraints has been presented 
recently in Ref. [10], including results from UA2 [31], the Tevatron [32], ATLAS [33, 34] and 
CMS [35]. These searches can set bounds on a broad range of masses and are complementary 
to the signatures we have studied in this paper. We will include the di-jet limits in the final 
discussion of the bounds on the model parameter space presented in section 5. 

4 Astrophysical and cosmological constraints 

Models for dark matter are constrained by the relic density, indirect and direct searches. We have 
evaluated the relic density and the most recent IceCube limits on dark matter annihilation in 
the Sun [23] within our simplified model for Majorana dark matter, as discussed in detail below. 
We also briefly comment on direct dark matter searches, which provide further complementary 
constraints on our model. 

4.1 Relic density 

We have calculated the dark matter relic density for each point of the model parameter space 
using FeynRules 2.1 [27], CalcHEP [36] and micrOMEGAs 4 [37]. The calculation includes 
possible resonant effects, which are not taken into account when expanding in the dark matter 
velocity [38]. We have however checked that the results obtained from the exact calculation 
within the simplified model match those obtained in the EFT limit for My m-^. To that end, 
we have calculated the annihilation cross section into quarks [39, 40] for the case of Majorana 
dark matter. 


<yqq — 


9g9l 1 /3g 


^4My [28mgm^ — 4s (m^ -|- m^) -|- s^] 
— 96Mym^mgS -|- 48m^mgS^^ , 


( 3 ) 


where = y 1 ~ and s is the center-of-mass energy. Note that this annihilation cross 

section is helicity suppressed [41, 42], and thus the main annihilation channels will be into top 
and bottom quarks. 

In Fig. 6 (grey band) we show the parameter values of our model for which the relic den¬ 
sity from thermal freeze-out agrees with the one measured by the Planck Collaboration [43], 
= 0.1199, within ±10%. Assuming a standard cosmological history parameter points 
above this line can be considered excluded in the framework of our simplified model. Points 
below this line could be allowed by either requiring an additional component of dark matter or 
an additional (non-thermal) production mechanism. In this paper we assume 100% of the (local) 
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dark matter to be the considered WIMP candidate and hence do not rescale the limits from 
IceCube and LUX for points below the grey band. Note that in an extension of our simplified 
model co-annihilations or additional resonances may exist that could further weaken the relic 
density constraint while providing a similar phenomenology at the LHC and IceCube. 

Finally, we mention that the unitarity of the S'-matrix imposes constraints on the masses and 
couplings in the dark matter model [10, 44]. However, performing a definite analysis of unitarity 
constraints would require us to extend the simplified model framework, which is beyond the 
scope of this work. 


4.2 Limits from IceCube 


If dark matter particles scatter in heavy astrophysical objects such as the Sun, they can lose 
enough energy to become gravitationally trapped inside the object. With the accumulation of 
dark matter, the annihilation rate can become large enough to lead to an equilibrium between 
dark matter capture and annihilation. 

The evolution of the number of dark matter particles in the Sun, N, can be described by the 
Riccati differential equation [45] 

N = Cq- CaN^ - CeN, (4) 


where N denotes the time derivative of N, Cq is the capture rate of dark matter particles in the 
Sun, CaN'^ = 2^4 is twice the dark matter annihilation rate and CeN is the evaporation rate, 
i.e. the rate at which particles escape the Sun due to hard elastic scattering. For dark matter 
particles with masses > lOGeV, the evaporation term can be neglected [46], allowing for a 
simple solution of Eq. (4): 

CaN^ = Cq tanh^ ( t) . (5) 

For large times, \JCqCa 1, the tanh-term in Eq. (5) approaches one, and WIMP annihilation 

and capture are in equilibrium, i.e. Cq = 2 T 4 = CaN'^- This implies = 0. Hence in 
equilibrium the annihilation rate does not depend on the annihilation cross section, but only on 
the capture rate, which in turn is determined by the elastic WIMP scattering cross section. We 
shall analyze the equilibrium condition within our model in more detail in Sec. 4.2.3. 

Through a measurement of the neutrino flux, neutrino telescopes are sensitive to the WIMP 
annihilation rate. We consider data from the IceCube Neutrino Observatory [47, 23] taken 
during 317 days in the years 2011 and 2012. No significant excess over background has been 
observed and these measurements can thus be used to set limits on possible dark matter signals. 

The search has been interpreted in terms of limits on the spin-dependent WIMP-proton 
scattering cross section. The WIMP-nucleon scattering cross section in our model is 
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(N) 

SD 


TVMy 




~ 1.8 X 10"'^°cm^ 



( 6 ) 


where ^aix — /i'^x + reduced WIMP-nucleon mass and N = p for WIMP- 

proton scattering.^ In the second line of Eq. (6) we have used the fact that pq is universal 
and included the numerical values for the nucleon form factors = A^^'^ = 0.85, A^'^ = 

^The scattering cross section for Majorana fermion dark matter given in Eq. (6) is larger than the one for 
Dirac fermion dark matter by a factor of four (see e.g. Ref. [48]). This is in contradiction with the result quoted 
in Ref. [49], which is a factor of eight smaller than our result. 
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= —0.42 and = A^'^ = —0.08 [50]. Considering universal couplings and neglecting the 
small mass difference between the proton and the neutron, the WIMP-neutron scattering cross 
section probed in direct detection experiments (see section 4.3) is equal to the WIMP-proton 
cross section. 

The dark matter interpretation of the IceCube searches initially relied on two scenarios: an¬ 
nihilation into bb and W~^W~ [47]. Recently, an improved interpretation [23] has been performed 
which includes more scenarios, in particular annihilation into tt. In the model considered here, 
annihilation into tt, bb and VV is dominant (c/. Fig. 6). Therefore, in Sec. 4.2.1 we will first 
estimate the limit on the spin-dependent WIMP-proton scattering cross section for annihilation 
into tt based on the limits for the annihilation into W~^W~ and bb. The excellent agreement be¬ 
tween our estimate of the tt limits and the results presented in Ref. [23] is considered as evidence 
for the accuracy of the conversion method. We then apply the same method to estimate limits 
on the annihilation into VV in Sec. 4.2.2. For the exclusion limits on the parameter space of 
our model we conservatively take into account only the most constraining channel and compute 
the model rejection factor, /r, via 


My) = My) x maxj 


Rti{mx,My) Ri,i{mx,My) Ryy{mx,My) 


af^irrix) 


O'; 


UL 

bb 


{rux) ’ a]^\r{mx,My) 


(7) 


where Ri is the contribution of the channel i to the annihilation cross section and is the 
corresponding upper limit on the spin-dependent WIMP-proton scattering cross section. Note 
that we display the dependence of the various factors on the dark matter and mediator masses, 
but not the dependence on the couplings gq,gx- A point in parameter space is excluded if /i > 1. 


4.2.1 Limit on cjsd for annihilation into tt 

In this subsection we will derive a limit on the spin-dependent WIMP-proton scattering cross 
section for dark matter annihilating to 100% into ti. Our procedure uses the limits of capturing 
rates 1% for annihilation into W^W~ and bb as input. As WIMPs are practically at rest inside 
the Sun, the W bosons (or b quarks) that are produced directly in WIMP annihilation have 
a well-defined energy given by the WIMP mass nix- Therefore we can interpret the limits 
on /a as limits on the annihilation rate of the IF-(or 6-)pairs with the energy 
regardless of the actual annihilation process.^ We exploit this fact and consider the energy 
spectrum of the W bosons and b quarks arising from the decay of the top quarks from WIMP 
annihilation into tt. Since the interaction time of the top quark in the Sun is much larger than 
its lifetime, no energy loss is expected before it decays into an on-shell VF-boson and 6-quark. 
As a next step we calculate the probability distributions P{Ei\Et) of the energy of a final state 
particle i for a given top energy - and hence WIMP mass - by simulating the annihilation 
process yy —)• tt —)■ W^W~bb with Madgraph5_aMC@NLO. The probability distributions 
are normalized to one. Based on these distributions we can calculate the resulting limit on 
by a weighted average of the limits for each relevant final state x = , W~, 6, 6: 


1 

niEt) 



P{Ew\Et) 

rr^-{Ew) 


+ 



P{Eb\Et) 

pfm 


( 8 ) 


In order to evaluate E^ ^ for arbitrary values of the energy we interpolate the limits 
linearly in E^^^ between the values given in Ref. [23] on a double logarithmic scale. Note that 

^Note that this is only true for the limit on Fa, but not for the limit on ctsd because ctsd depends explicitly 
on through the capture efficiencies. Furthermore, the limit on Fa is independent of the correlation between 
the two W bosons (or b quarks) of one annihilation. The fraction of events where multiple neutrinos arising from 
W bosons (or b quarks) of the same annihilation simultaneously interact in the detector can safely be neglected. 









the contribution from the b quarks in Eq. (8) is sub-leading. The correction is below 10% (for 
large WIMP masses and decreases to below 1% towards small masses). 



Figure 4: 90% CL upper limits on the spin-dependent WIMP-proton scattering cross section, ctsdi lor 100% 
annihilation into tt (red solid line), bb (blue dotted line) and W^W~ (green dashed line) from the IceCube 
Neutrino Observatory [23]. The black dashed line shows the tt limit we derived from the bb and W'^W~ data 
according to the method described in Sec. 4.2.1. 

We converted the limits on the spin-dependent WIMP-proton scattering cross section given 
in Ref. [23] into limits on the annihilation rate (and vice versa) by 

= ^^SD(m^) (?7tx), (9) 

where Rsd is the (channel independent) capture efficiency for the spin-dependent part of the 
scattering. We take Rsd from Ref. [47] where it was derived following the method of Ref. [51] 
adopting the Standard Solar Model BS05(OP) [52] as implemented in DarkSUSY [53]. 

The result is presented in Fig. 4, where we also show the limit on 100% annihilation into 
W~^W~, bb and tt from Ref. [23]. The good agreement of our cross section limits for tt with 
Ref. [23] demonstrates the validity of our calculation. In addition, as a consistency check of 
our procedure, we have converted the limit on the annihilation rate into a limit on the resulting 
muon flux in the detector, using the conversion functions provided in Ref. [51]. As expected, 
we found that from the annihilation into tt is always in between the limits obtained for 
annihilation into W^W~ and bb as presented in Ref. [23]. 

4.2.2 Limit on cjsd for annihilation into VV 

In this subsection we will derive limits on the spin-dependent WIMP-proton scattering cross 
section for dark matter annihilating to 100% into mediator pairs VV with the method described 
in Sec. 4.2.2. Annihilation into VV takes place only for My < m^, where the mediator decays 
solely into quarks. As we consider a universal coupling to all quarks, the corresponding branching 
ratios are simply determined by the accessible phase space. Here we only take into account 
neutrinos arising from the decay of the mediator into bottom and top quarks. In order to 
justify this choice we computed the differential neutrino spectra driu/dEi^ for annihilation into 
all quark flavors with WiMPSiM 3 [54].“^ As expected, the neutrino fluxes spectra for light flavor 
quarks d, u, s are much softer than for bottom and top quarks and can be safely neglected for 
the derivation of the limits. The neutrino flux for annihilation into charm quarks is weaker than 

^The program package WimpSim is linked to PYTHIA 6 [29] for the simulation of dark matter annihilation 
in the Sun, Nusigma [55] for the simulation of neutrino-nucleon interactions and to DarkSUSY [53] for the 
implementation of the Sun’s density profile. We take the neutrino oscillation parameters from Ref. [56]. 
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that for 6-quarks by a factor of 3 to 10 (and much weaker than for top quarks) in the relevant 
energy range and is hence subdominant. We simulate the annihilation process yx arid 

the subsequent decays V ^ qq with Madgraph5_aMC@NLO and determine the probability 
distributions P{Eq\Ev) of the energy of a final state quark q = b,t, which are hence normalized 
to 2 X BR(y — )• qq). We calculate the resulting limit on E^^ analogous to Eq. (8) through: 


1 

rX^iEv) 



P{Et\Ev) 

EXiEt) 


+ 



P{Eb\Ev) 

PXiEb) 


( 10 ) 


As the resulting limit depends on the WIMP mass and the mediator mass we scan the 
corresponding two-dimensional grid. In Fig. 5 we present the limits on Ia (left panel) and 
(right panel). The limits are considerably weaker than for annihilation into a pair of tops. As 
the limits depend both on rriy’ and My we show two slices in the parameter space My/m^ = 0.75 
and My/rriy^ = 0.35, respectively. The main difference between these two slices is due to the 
opening of the mediator decay into top quarks for My > 2mt, which greatly enhances the 
sensitivity. 


The resulting IceCube limits on the model parameter space considering the annihilation 
channels 66, ti and VV are shown in Fig. 6. The sensitivity to annihilation into VV is signifi¬ 
cantly weaker than for tt and - below the tt threshold - for 66 final states. This causes a drop 
in the limit on My for regions where annihilation into VV is dominant (light grey areas). 




Figure 5: 90% CL upper limits on the annihilation rate Ta and the spin-dependent WIMP-proton scattering 
cross section ctsd for 100% annihilation into VV (black lines) derived from a reinterpretation of the data from the 
IceCube Neutrino Observatory [23]. We show two slices regarding the mediator mass, Mv/m^ = 0.75 (dashed 
line) and My jm^ = 0.35 (dot-dashed line). For comparison we also show the limits for 100% annihilation into bb 
(blue solid line) and W'^W~ (red solid line) taken from Ref. [23]. 


4.2.3 Equilibrium condition for capturing and annihilation in the Sun 


As discussed at the beginning of Sec. 4.2, for large times, ^/CQCA t !§> 1, dark matter matter 
annihilation and capture in the Sun are in equilibrium, and the limits on the annihilation rate 
can directly be translated into limits on the elastic WIMP-nucleon scattering. Assuming that 
the Sun has been collecting dark matter during its entire lifetime, t = tQ ~ 1.5 X 10^’^ s, the 
equilibrium condition can approximately be expressed by [45] 


^C^CaIq ~ 330 



{aAv) 

cm^s"^/ VlOGeV 


3/4 

> 1 . 


( 11 ) 


In practice \JCqCa t© ^ 3 is already enough to obtain an error of less than a percent on C©. 
In order to estimate ^JCqCa t© for our model, we compute (oav) with micrOMEGAs. The 
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Figure 6: Astrophysical and cosmological quantities in the m^-My plane in four slices of the considered parameter 
space: = 1, Fy = O.OlMy (upper left panel), yjQxQq = 1, Fy = 0.5My (upper right panel), x/gxQq = 0.2, 

Fy = O.OlMy (lower left panel) and ^/Qxdq = 0.2, Fy = 0.5My (lower right panel). The shaded regions denote 
the dominant annihilation channel, in red, blue, green and grey we denote dominant annihilation into ti, bb, light 
flavor quarks and two mediators, respectively. The 90% CL lower exclusion limits from the IceCube Neutrino 
Observatory (green line) and from LUX derived from the limits presented in Ref. [49] (purple line) are also 
displayed. The dark grey shaded band denotes the region where the relic density matches the dark matter density 
within ±10%. In the region below the red thin line the equilibrium condition is fulfilled, i.e. \/CqCa tQ > 3. 


capture rate, Cq, can be deduced from the WIMP-nucleon scattering cross section using the 
capture efficiency (see Sec. 4.2.1) 


Cq = KsB{'rny^)crsD ■ ( 12 ) 

In Fig. 6 we show the contours x/CqCa tQ = 3 (labeled “EQ”) in the m^-My plane. Below this 
line the equilibrium condition is fulfilled and the interpretation of the IceCube measurement in 
terms of the elastic scattering cross section is justified. A similar conclusion was put forward 
for the limits set in Ref. [4], where {crjxv) was computed within the EFT. However, the fact 
that the equilibrium condition holds in the EFT limit does not imply that it should hold in the 
simplified model description, as the thermally averaged cross section, (uau), can be smaller in 
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the latter case. 


4.3 Limits from direct detection 

Although providing much stronger limits on the spin-independent WIMP-nucleon scattering 
cross section, direct detection experiments are also capable of imposing limits on spin-dependent 
scattering. In the relevant part of the parameter space the strongest constraints are set by the 
LUX experiment [57]. As an up-to-date dedicated analysis for spin-dependent scattering has 
not been provided by the LUX Collaboration, several authors have reinterpreted the LUX limits 
accordingly [7, 49, 10]. 

We will use the limits on the spin-dependent WIMP-neutron scattering cross section pre¬ 
sented in Ref. [49]. They are very similar to those presented in Ref. [7]. Since Xenon, the target 
material of LUX, has neutron-odd isotopes, limits on spin-dependent WIMP-proton scattering 
are considerably weaker and are thus not taken into account. The results are shown in Fig. 6 
(purple curve). 

5 Results and Discussion 

In this section we finally summarize the constraints on our model from indirect, direct and 
collider searches. Fig. 7 presents limits on the dark matter and mediator masses for four different 
choices of the coupling strength, y/gx9q^ a-iid the mediator width, Fy. As discussed in Sec. 2, 
not all combinations of the parameters m^. My, Fy and yjQxdq ^re viable within our model. For 
small width, a large range of the parameter space, marked in shaded blue in the left panels of 
Fig. 7, is not allowed. For a large width, Fy = 0.5 My, on the other hand, the whole parameter 
region is allowed. 

The model provides the correct dark matter relic density only for the parameter region within 
the thin grey band. The shaded grey area in Fig. 7 indicates the parameter region leading to 
an overproduction of dark matter within the model. Note, however, that the relic density 
constraint can be softened significantly if we assume an extended particle spectrum that leads 
to co-annihilation effects (which can give access to the region above the thermal relic-band) or to 
additional - non-thermal - contributions to the dark matter production (which can give access 
to the region below the thermal relic-band). 

The limits from the ATLAS and CMS mono-jet searches, obtained within the simplified 
model, are particularly relevant for small rux- The ATLAS mono-jet analysis (blue solid curve) 
is the most constraining search for My > 2m^ in the whole parameter space considered. In 
particular, this search constrains most strongly the parameter space where the relic density 
from thermal freeze-out agrees with the measured dark matter density, c/. the dark grey band 
that represents the parameter region with = 0.1199 ±0.012. The LHC limits depend upon 

the width of the mediator. This dependence is particularly pronounced for -yjQxSq = 0-2, where 
LHC searches probe mediator masses My < 1 TeV. Note that this mass range is of the order of 
the typical scattering energies of mono-jet searches with large MET, so that an EFT description 
would not be reliable. The sensitivity of the LHC searches is higher for smaller width as the 
cross section is enhanced through on-shell mediator production. 

As discussed in Sec. 3.2 we also show the limits from di-jet searches as obtained in Ref. [10]. 
These constraints are particularly important for large mediator width (right panels), because a 
larger width requires a larger mediator-quark coupling gq^ and thus a larger cross section. 

The direct and indirect dark matter searches by the LUX and IceCube experiments, respec¬ 
tively, are particularly relevant for intermediate and large dark matter masses. The IceCube 
limits, specifically, show a maximal sensitivity for dark matter masses around 200—1000 GeV. 
This results from two effects: As the sensitivity for annihilation into bb is significantly lower 
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[GeV] 


x/9x9q = 0-2, rv = O.OlMy 



x/9x9q = 0-2, rv = 0.5Mv 



Figure 7: Summary of the exclusion limits in the m^-Mv plane in four slices of the considered parameter space: 
y/gxQq = 1, Fy = O.OlMy (upper left panel), y/fhdh = 1, Fy = 0.5My (upper right panel), y/g^Qq = 0.2, 
Fy = O.OlMy (lower left panel) and = 0.2, Fy = 0.5My (lower right panel). We show the 95% CL 

lower exclusion limits from mono-jet searches at ATLAS (blue lines) and CMS (red lines) (both in the simplified 
model interpretation) as well as limits from searches for resonances in di-jet signatures taken from Ref. [10] (orange 
shaded regions are excluded). Furthermore, we show the 90% CL lower exclusion limits from the IceCube Neutrino 
Observatory (green line) and from LUX derived from the limits presented in Ref. [49] (purple line). The dark 
grey shaded band denotes the region where the relic density matches the dark matter density within ±10%. In 
the light-grey shaded region above it, the dark matter is over-produced. The blue shaded region in the left panels 
do not allow for a consistent solution for the mediator width as a function of Mv,m^, ^gxgq within the model. 


than for tt {cf. Fig. 4) there is a significant strengthening of the IceCube limits around the top 
threshold. For very large on the other hand, the capture efficiencies and thus the cross 
section limits decrease and the IceCube search loses sensitivity. Furthermore, for My < 
the annihilation into a pair of mediators is allowed (via a t-channel x), which again provides a 
lower model rejection than tt. Hence, the limit depends strongly on the relative importance of 
annihilation into tt and VV. For > My and x/g^Jjq = 1, the limit is considerably weaker 
for a small mediator width, as a smaller width requires a smaller gq and hence a larger (note 
that the cross section for annihilation into mediators scales like g^). 
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6 Conclusion 


The complementarity of direct, indirect and collider searches for dark matter can be exploited 
in the framework of simplified models, where dark matter and its experimental signatures are 
described with a minimal amount of new particles, interactions and model parameters. We 
have considered a simplified model with a Majorana fermion dark matter particle and an axial- 
vector mediator with universal couplings to SM quarks. Such a model leads to spin-dependent 
interactions and can thus be probed by IceCube in the search for neutrinos from dark matter 
annihilation in the Sun. 

We have focused on a re-interpretation of IceCube limits on the dark matter annihilation rate 
within our model, and have obtained substantial constraints on the dark matter and mediator 
masses. Furthermore, we have derived new constraints on our model from recent LHC mono¬ 
jet searches, and have analyzed in some detail the differences between interpretations of LHC 
searches within simplified models and effective field theories. 

We find that the limits from the ATLAS and CMS mono-jet searches are particularly relevant 
for small dark matter masses. They exclude mediator masses My < 1 TeV, depending in detail 
on the size of the couplings and the mediator width. The indirect searches for dark matter 
annihilation in the Sun by IceCube probe intermediate and large dark matter masses and show 
a maximal sensitivity for masses ~ 200—1000 GeV. In this region, the dark matter capture 
efficiencies in the Sun are still sizeable, and dark matter annihilation is predominantly into top 
quarks, leading to more highly energetic neutrinos and thus a higher neutrino detection efficiency 
with IceCube. 

We have also computed the relic density within our model, and have combined the LHC 
mono-jet and IceCube limits with constraints from direct detection and the collider search 
for di-jet resonances. We have found a striking complementarity of the different experimental 
approaches, which probe particular and often distinct regions of the model parameter space. 
Thus, the combination of future collider, indirect and direct searches for dark matter will allow 
a comprehensive test of minimal dark matter models. 
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